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Abstract. In this paper, the necessary and sufficient conditions in 
order that a smooth mapping F(r, a, /3, a, b) be a dependence of a com- 
plete solution x(r) of some second-order ordinary differential equation 
on Neumann conditions x(a) = a, x(/3) — b, a 7^ (3 are deduced. These 
necessary and sufficient conditions consist of functional equations for 
F and of a smooth extensibility condition. Illustrative examples are 
presented to demonstrate this result. In these examples, the mentioned 
functional equations for F are related to the functional equations for 
geodesies, to Jensen's equation, to the functional equations for conic 
sections and to Neuman's result for linear ordinary differential equations. 
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1. Introduction 

The relationship between functional and ordinary differential equations was 
studied by many authors, e.g., M. A. Abdelkader [1], E. Castillo and R. Ruiz- 
Cobo [4], L. N. Eshukov [5]. We will discuss here a functional formulation 
of second-order ordinary differential equations. 

In Section 2, we present the functional equations (1), (2) for a continuous 
function of many variables. A solution for 1-dimensional linear case is found. 
This solution is related to Neuman's result for linear second-order ordinary 
differential equations [10]. 

In Section 3, we investigate the relationship between the introduced func- 
tional equations and ordinary differential equations in general. A notion of 
so-called N-solvable second-order ordinary differential equation is presented. 
We prove that under given condition (12) the functional equations are equiv- 
alent to the family of all smooth N-solvable second-order ordinary differential 
equations. 

Section 4 contains examples. The introduced functional equations are 
related to Jensen's functional equation, to the functional equations for 
geodesies, and also to the functional equations for conic sections. 

For the sake of simplicity, in this paper we consider smooth ordinary 
differential equations globally defined on the Cartesian product and globally 
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N-solvable. 



2. Functional equations 

Let n be a natural number, D = {(a,f3) € M 2 | a = (3} the diagonal of M 2 . 
Let F: M. x (M. 2 \ D) x M 2n — > M n be a continuous mapping. We consider 
functional equations 

F(t, a, (3, a, b) = F(r, 7 , 5, F( 7 , a, (3, a, b), F(S, a, f3, a, b)), (1) 

F(a, a, (3, a,b) = a, F((3,a, (3,a,b) = b. (2) 
Let us start with an example. 

Example 1 (linear case). Let n = 1 and let F be linear in the last two 
arguments. Then F is a solution of the functional equations (1), (2) if and 
only if there exist functions x\,X2- R — > K such that whenever a (3 the 
following statements hold 

xi(a)x 2 (/3) - xi(/3)x2(a) / 0, (3) 

F(r,a,(3,a,b) = 

(x 1 (t)x 2 ((3) - xi{(3)x 2 {T))a + (xi(a)x 2 (r) - xi(r)x 2 (a))^ 
xi (a)x2 {(3) - x 1 {(3)x2 (a) 

Let us check it. Conditions (1), (2) can be obtained from (3), (4) by direct 
computation. Conversely, put x\(t) = F(t,0, 1, 1,0), x 2 (r) = F(t, 0, 1,0, 1). 
Then 

(x\(a) x 2 (a) a \ 

xi (J3) x 2 (J3) b =0. (5) 

xi(r) x 2 (r) F(T,a,/3,a,b) J 

Indeed, from linearity and from (1) the third row is a linear combination 
of the first row and the second row with coefficients F(r, a, (3, 1, 0) and 
F(t, a, (3, 0, 1). Since from linearity and from (1), (2) 

( F(0,a,/3,1,0) F(0,a, (3,0,1) \ f Xl (a) x 2 (a) \ / 1 0\ 
^ F(l,«,/3,1,0) F(l,a,/3,0,1) A xi(/3) ^(/3) J v 1 J 

holds, we have the inequality (3). Finally, from (5) we can express F in the 
form (4). 

Remark 1. Consider a second-order linear differential equation with linearly 
independent solutions xi,x 2 . For this case the condition (5), where F is 
the solution satisfying (2), was presented (in a more general form) by 
F. Neuman [10] in 2003. In the following chapter we will investigate the 
relationship between the functional equations (1), (2) and differential 
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equations in general. 



3. The relationship between functional and differential 

equations 

Let /:lxR"x R n — > R ra be a mapping. Consider a second-order equation 

x = f(r,x,x) (6) 

and Neumann conditions 

x(a) = a, x((3) = b, (7) 

where a, (5 G R, a / /3, a, 6 G R n . 

The mapping F will be called a dependence of the solution on Neumann 
conditions if and only if the following equality holds 

F(T,a,(3,x(a),x(/3))=x(T) (8) 

for any a / /3, r G R and for any complete solution x of the equation 
(6). Notice that in this case Dom(x) = R. We say that the dependence of 
the solution on Neumann conditions is smooth if and only if F is a smooth 
mapping. The equation (6) is called N-solvable if and only if there exists 
a unique complete solution for any Neumann conditions (7). Consider the 
equation (6) with integral conditions 

l 

x(a) = a, J x((l — j)a + 7/3) d-f = v, (9) 
o 

where a,/3,j G R, a,v G R n . The equation (6) is called I-solvable if and 
only if there exists a unique complete solution for any integral conditions. 

Lemma 1. The smooth equation (6) is N-solvable if and only if it is 
I-solvable. 



Proof. Let (6) be N-solvable and smooth. For a ^ (3 we can rewrite (9) as 

x(a) = a, x({3) = a + v((3 — a). (10) 

For a = (3 we have (9) in the form 

x(a) = a, x(a) = v. (11) 

The existence and uniqueness of the solution of (6), (10) (respectively, (6), 
(11)) follows from the assumption of N-solvability (respectively, from the 
uniqueness of the solution of the Cauchy problem, see, e.g., [6]). Therefore 
(6) is I-solvable. 

Conversely, let (6) be I-solvable. Then (6) is N-solvable. This follows 
from (10) immediately. □ 
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Lemma 2. The dependence (r, a, j3, a, v) i— > x(r) of the solution x of the 
smooth I-solvable equation (6) on conditions (9) is smooth on an open set 
containing txDx R 2n , where D is the diagonal of R 2 . 

Proof. Let r i— > x(t) be the solution of (6) satisfying Cauchy conditions 
x(a) = a, x(a) = u. The dependence of x on r,a,a,u is smooth (see, e.g., 
[6, chapter V]). The mappings 



are smooth too. It remains to check that the dependence of u on a, (3, a, v is 
smooth. If we put @ = a in (9), using Cauchy conditions we get u = v. Since 
h(oto, ao, ao, vq, vq) = and d(ag, oq, ao, vq, vq) = 1 for ao € R, ao,t»o £ R n , 
we can apply [9, Theorem 1.3.5. and Corollary 1.3.9.] to the mapping h. 
From this there exists a neighbourhood U\ C R 2n+2 of (ao, ao, ao, vq) and 
a neighbourhood U2 C M n of «o such that for any (a,f3,a,v) € U± there 
is a unique u = u(a, (3,a,v) € U2, which satisfies h(a, (3,a,v,u) = and 
(a, /?, a, v) 1— > u(a, /?, a, t> ) is a smooth mapping. Finally, the dependence of 
x on r and on parameters (a, /?, a, v) from E/i is smooth. □ 

Theorem. A mapping F is the smooth dependence of the solution of some 
smooth N-solvable equation (6) on Neumann conditions if and only if the 
conditions (1), (2) hold and a mapping 



has a smooth extension S to R 2n+3 . 

Proof of Theorem. Assume that F is the smooth dependence of the solution 
of some smooth N-solvable equation (6) on Neumann conditions (7). From 



(7), (8) we get F(r, a, (3, a,b) = x(t) = F(T,^,5,x( / y),x(5)) for a 7^ (3, 



7 7^ 5, which implies (1). If we put a = r (respectively, f3 = r) in (8), using 
(7) we obtain first equality in (2) (respectively, second equality in (2)). We 
define a mapping S: (r, a, (3,a,v) 1— > x(t), where x is the solution of the 
smooth N-solvable equation (6) satisfying conditions (9). The mapping S 
exists for any (r, a, f3,a,v) € M 2n+3 , this follows from Lemma 1. Moreover, 
from (8), (10) we have S(t, a, (3, a, v) = F(t, a, (3,a,a + v((3 — a)) for a 7^ (3. 
Therefore S is the extension of (12) to IR 2n+3 . The smoothness of S follows 
from the smoothness of F and from Lemma 2. 



h: (a, (3, a, v, u) 1— > / x((l — 7)0 + 7/?) 07 — v, 








R x (R 2 \D) x R 2n 3 (r, a, (3, a, v) ^ 

F(t, a, (3,a,a + v{(3 - a)) G R' 



(12) 
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Now we assume that (1), (2), (12) hold. With help of (12) we calculate 
the limit 5 —>■ 7 in (1) as 

dS^, a, (3, a, v) ' 



S(t, a, P,a,v) = S[ r, 7, 7, S{j, a, (3, a, v), 
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Differentiating this with respect to r two times and putting (3 = a, 7 = r 
we see that r 1— > 5(r, a, a, a, v) is the solution of (6), where 

d 2 S(r, 7, 7, x, x) 



f(r,x,x) 



8t 2 



(13) 



7=T 



From (2), (12) we obtain S(a,a,a,a,v) = a and 



8S(t, a, a, a, v) 



dr 



8S(t, a, (3, a, v) 



dr 



j3=a 



d_ 

dp 

d_ 

~~ 3(3 



(s(T,a,/3,a,v) \ T= p) 



dS(r, a, (3, a, v) 



dp 



T=0j 



(3=a 



(s(T,a,/3,a,v)\ T=/3 - S(r, a, (3, a, v)| T=Q ) 



(3=a 



— (a + v(P - a) - a) 



v. 



From (13) / is smooth. Let x be a complete solution of (6), (13). From 
above we can see that the solutions r 1— > S(t, a, a, x(a),x(a)) and r 1— > x(t) 
satisfy the same Cauchy condition. Therefore 



S(T,a,a,x(a),x(a)) = x(t) 



(14) 



for any t£R (see, e.g., [6, chapter V]). Finally, as (3 — > a in (1), using (7), 
(12), (14) we obtain F(t,j,6,x('j),x(8)) = x(t), which is equivalent to (8). 
Hence F is the dependence of the solution on Neumann conditions. The 
smoothness of F follows from (12). The uniqueness of the complete solution 
of (6), (7), (13) follows from (8). Thus the equation (6), (13) is N-solvable. □ 



Remark 2. From Theorem we can see that under the condition (12) the func- 
tional equations (1), (2) are equivalent to the family of all smooth N-solvable 
ordinary differential equations (6). It follows from the proof of Lemma 2 
that every smooth differential equation (6) is locally I-solvable, i.e., for any 
ao € K, ao,uo £ W 1 there exists a neighbourhood Uq C K 2n+1 of (ao; a tb^o) 
and a neighbourhood U\ C M. 2n+2 of (ao, ao, ao> ^o) such that the equation 
(6) restricted on Uq has a unique complete solution for any (a, (3, a, v) £ U\. 
Moreover, from the proof of Lemma 1, every smooth equation (6) is locally 
N-solvable. Thus the functional equations (1), (2) hold locally for any 
smooth equation (6), but Dom(F) is more complicated in such a case. 
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4. Examples 

Example 2 (free fall). Let us consider the free fall equation x = g. It is 
easy to see, that for this equation we have 

F(r, a, (J, a, b) = ^ - P)a+(a - r)b + ^ _ ^ _ p) (lg) 

and a smooth extension (12) has the form 

(r, a, (3,a,v) i-> a + f (r — a) + ^(r — a)(r — /?). 



Example 3 (geodesies). A geodesic can be viewed as a solution of a second- 
order differential equation. Thus this is the special case of the general situ- 
ation presented above. Let T be a smooth linear connection on a manifold 
M. Let U C M be an open convex set and let a,b € U. By G ab we denote 
the geodesic of T satisfying G ab (0) = a, G ab (l) = b, G ab ([0, 1]) C U. The 
existence of U and uniqueness of G ab is guaranteed by the Whitehead local 
convex lemma (see, e.g., [8] or [11]). We define a mapping G in the following 
way 

G: (a,b,p) i-> G ab (p). 

Since there exist constant geodesies, we have 

G(a, a, p) = a. (16) 

Let x be a geodesic satisfying x(a) = a, x{(3) = b. Since every affine 
reparametrization of a geodesic is a geodesic again, we have x(t) = 
G(a, b, tjE^)- If there exists a dependence of the solution on Neumann con- 
ditions, then from the previous equality and from (8) we obtain 

F(r,a,P,a,b) =G^a,b,^^j , (17) 

where Dom(G) = l n xl n xt and Codom(G) = R n . The following functional 
equations for geodesies were deduced by L. Klapka [7] in 2000 

G(a,b,0)=a, G(a,b,l) = b, (18) 

G(a, b, (1 - p)( + m) = G(G(a, b, Q, G(a, b, v ), p). (19) 

Under the conditions (16), (17), the statements (18), (19) follow from (1), 
(2) for any (, rf, p G R. 

Example 4 (Jensen's functional equation). Consider a mapping G from 
Example 3. Let G be linear in the first two arguments. By putting £ = 1, 
rf = in (18), (19) we obtain G(a, b, 1 — p) = G(b, a, p). For any p G K we 
define a mapping Q(p) in the following way 



Q(p): R n Bb^G{0,b,p) £ R n . 
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From this and from linearity we get G(a,b,p) = Q(l — p){a) + Q(p)(b) 
and using (16) we obtain Q (j^J (a) = \a. Finally, by substituting a = 0, 

p = \ into (19) and using the previous equalities we get Jensen's functional 
equation 

Q fC + v\_ Q(C) + Q(v) 



Example 5 (Conic sections). Let n = 1. Let F be continuous, linear inho- 
mogeneous in the last two arguments, autonomous, i.e., 

F(t, a, (3, a, b) = F(t + A, a + A, (3 + A, a, b) (20) 

for any A G R and nondissipative, i.e., F(r, a, (3, a, b) = F(— r, —a, —(3, a, b). 
Under these assumptions we can find a solution of (1), (2). Let a,(3,a,b be 
fixed. We denote x(t) = F(r,a, (3,a,b). From (1), (20) we get 

x( T ) = F(0, A, 2A, x(t + A), x(t + 2A)), 

x{t + 3A) = F(0, -A, -2A, x{t + 2A), x(t + A)) 

for any r G R, A G R\{0}. By using the inhomogeneous linearity and 
the nondissipativity we see that the difference x(t + 3A) — x(r) depends on 
x(r + 2A) — x(r + A) linearly and a coefficient depends on A only. Similarly, 
x(t + 4A) — x(t + A) depends on x(t + 3A) — x(r + 2A) linearly and the 
coefficient is the same. Therefore 

(x(t + 4A) - x(t + A))(z(t + 2A) - x(t + A)) = 

(x(r + 3A) - x(t)){x(t + 3A) - x(t + 2A)). 

This is the functional equation for characterization of conic sections pre- 
sented by A. Angelesco (see [2, section 2.5.2.] or [3]) in 1922. All its continu- 
ous solutions are known and these solutions are real analytic. Differentiating 
this equation with respect to A five times and using the limit A^Owe get 
an equation 

zya rv* ry* rv* 

•J,' ib Ju ^ 

which is equivalent to 

— rm) = o. 

dr V x(t) ) 

From this x = Kx and by integrating we obtain 'x = Kx+g, where K, g £ R. 
According to N-solvability K > 0. Putting K = k 2 we get 

f{r,x,x) = k 2 x + g. 

From (6), (8) 

F(t, a, (3, a, b) = 

(k 2 a + g)smh(k(T - (3)) + (k 2 b + £r)sinh(A;(a - r)) - g sinh(A:(a - (3)) 

k 2 sinh(/c(a — (3)) 
where k, g G R, k ^ 0. For fc = 0we get the free fall (15). 
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